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Abstract 

We report on a search for the decays J/ip -^ Dje^Ve + c.c., J /'ijj -^ D~e^i'e + c.c., 
and J/^ -^ D'^e'^C + c.c. in a sample of 5.8 x lO'^ J/^p events collected with the 
BESII detector at the BEPC. No excess of signal above background is observed, 
and 90% confidence level upper limits on the branching fractions are set: B{J/iIj -^ 
Dje+Ve + c.c.) < 4.8xl0-^ S(J/V' -^ D- e+ Ue + c.c.) < 1.2xl0-^ and S(J/V' ^ 
L>Oe+e- + c.c.) < 1.1x10'^ 

PACS: 13.25.Gv, 13.30.Ce 



1 INTRODUCTION 



Hadronic, electromagnetic, and radiative decays of the J/ip have been widely 
studied. However there have been few searches for rare weak J/ip decay pro- 
cesses. Kinematically, the J/ip cannot decay into a pair of charmed D mesons, 
but can decay into a single D meson. Searches for weak decays of J/ip to single 
D or Ds mesons provide tests of standard model (SM) theory and serve as 
a probe of new physics [1], such as TopColor models, the minimal supersym- 
metric standard model with or without R-parity, and the two Higgs doublet 
model [2]. 

The branching fractions of J/ip decays to single D or Dg mesons are pre- 
dicted to be about 10~^ or smaller [3] in the SM. The flavor changing neu- 
tral current (FCNC) process c — *> m occurs in the standard model only at 
the loop level where it is suppressed by the GIM mechanism. Fig. 1 shows 
the dominant Feynman diagrams within the standard model for the decays 



J/ip -^ D~e^Ve-, J/'4' ~^ D~e^Vei and J/ip — > D^e^e~ . No decay of this 
type has been observed before. In this paper, we perform a search for the 
decays J /ip -^ D~e^Ve-, J/'4' ~^ D~e~^i'e, and J/ip -^ D^e^e~ in a sample 
of 5.8x10^7/-?/' events collected with the Beijing Spectrometer (BESII) [4] de- 
tector at the Beijing Electron- Positron Collider (BEPC) [5]. Throughout this 
paper the charge conjugate states are implicitly included. 
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Fig. 1. Typical Feynman diagrams for (a) J /^tp -^ Dg e^i'e, (b) J/ip ^ D e^Ve-, 
and (c) J/V'^^e+e". 



2 BESII Detector 



BES is a conventional solenoidal magnetic detector that is described in detail 
in Ref. [6]. BESII is the upgraded version of the BES detector [4]. A 12-layer 
Vertex Chamber (VC) surrounds the beryllium beam pipe and provides track 
and trigger information. A forty-layer main drift chamber (MDC) located just 
outside the VC provides measurements of charged particle trajectories over 
85% of the total solid angle; it also provides ionization energy loss (dE/dx) 
measurements which are used for particle identification (PID). A momentum 
resolution of 1.78% VI + p'^ {p in GeV/c) and a dE/dx resolution for Bhabha 
electrons of ~8% are obtained. An array of 48 scintillation counters surround- 
ing the MDC measures the time of fiight (TOF) of charged particles with a 
resolution of about 200 ps for hadrons. Outside the TOF counters, a 12 radia- 
tion length, lead-gas barrel shower counter (BSC), operating in self quenching 
streamer mode, measures the energies and positions of electrons and photons 
over 80% of the total sohd angle with resolutions of cte/E = 0.21/^/E {E in 
GeV), 0-0 = 7.9 mrad, and (7^ = 2.3 cm. Outside the solenoidal coil, which 
provides a 0.4 T magnetic field over the tracking volume, is an iron fiux re- 
turn that is instrumented with three double-layer muon counters that identify 
muons with momentum greater than 500 MeV/c. 



In this analysis, a GEANT3-based Monte Carlo program (SIMBES) with de- 
tailed consideration of the detector performance is used. The consistency be- 
tween data and Monte Carlo has been checked in many physics channels from 
J /ip and ip{2S) decays, and the agreement is reasonable, as described in detail 
in Ref. [7]. 



3 Event Selection 



Each charged track is required to be well fitted to a helix that is within the 
polar angle region |cos^| < 0.8 and to originate from the beam interaction 
region, which is defined to be within 2 cm of the beam line in the transverse 
plane and within 20 cm of the interaction point along the beam direction. For 
particle identification, confidence levels {CL) are calculated for each particle 
hypothesis using combined time-of-fiight [8] and MDC energy loss information. 
Pions and kaons are identified by requiring the confidence level for the desired 
hypothesis to be greater than 0.1% and, further, by requiring the normalized 
weights, defined as C L^/ {C L^^+C Lk) , where a denotes the desired particle, to 
exceed 0.7. For electron identification, the dE/dx, TOF, and BSC information 
are combined to form a particle identification confidence level for the electron 
hypothesis. An electron candidate is required to have CLe > 1.0% and satisfy 
CLJiCLe + CL^ + CLk) > 0.85. 

An isolated neutral cluster is considered to be a photon candidate when the 
angle between the nearest charged track and the cluster is greater than 18°, 
the first hit is in the beginning six radiation lengths, the difference between the 
angle of the cluster development direction in the BSC and the photon emission 
direction is less than 37°, and the energy deposit in the shower counter is 
greater than 60 MeV. 



3.1 J/ip —^ D'e+Ue and J/ip —^ ^"e+z/e 



Both J/ip -^ D~e^Ve and J/ip -^ D~e^Ve candidate decays must contain 
an e"*". D~ mesons are reconstructed in two modes D~ -^ (/)7r~ and K~K*^, 
with (f) -^ K^K~ and K*^ — » K^ti' . The candidates are reconstructed 
from two oppositely-charged kaons and must have an invariant mass \Mxk ~ 
1.02| < 0.015 GeV/c^. The K*^ candidates are constructed from K^ and vr" 
candidates and are required to have an invariant mass in the range \Mk-k — 
0.8961 < 0.060 GeV/c2. 

D^ candidates are reconstructed in the mode D^ — ^ K^-k^ti^ . Candidate 
events are required to have four charged tracks which satisfy charged track 



selection criteria, and the total charge of the tracks is required to be zero. 

The polarization of the K*^ meson in Dj decay is also utilized to reject back- 
grounds by imposing a selection requirement on the helicity angle Oh- The 
helicity angle is defined as the angle between one of the decay products of 
the K*^ and the direction of the flight of K*'^ in the K*^ rest frame. Back- 
ground events are distributed uniformly in cos 6h since they originate from 
random combinations, while signal events are distributed as cos^Oh- The K*^ 
candidates are required to have Icos^j^l > 0.4. Fig. 2 shows the invariant 
mass distributions of K^K~ and K^ti^ systems with arrows indicating the 
selection of (p and i^*(892) candidates. 
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Fig. 2. (a) K^K~ and (b) K^tt invariant mass distributions of 
D- -^ (pTT-, (t> ^ K+K- and D' -^ K'K*^, K*^ -^ K+tt' candidates. 
The signal regions for individual channels are indicated by arrows in the plot. 

The angle between the identified electron and the nearest charged track is 
required to be larger than 12° to reject backgrounds from gamma conver- 
sions and 7r° Dalitz decays. In order to reject backgrounds from J/ip decaying 
to states with extra neutral particles, the number of isolated photons is con- 
strained to be zero. The requirements for an isolated photon are more stringent 
than those for a good photon in order to increase the detection efficiency. An 
isolated photon is a photon with the angle between the nearest charged track 
and the cluster greater than 22° and the difference between the angle of the 
cluster development direction in the BSC and the photon emission direction 
less than 60°, and the energy deposit in the shower counter is greater than 
0.1 GeV. 



The four-momentum of the neutrino, 



is inferred from 



Pu y^miss) Pmissjy 

the difference between the net four momentum of the J/ip particle, pj/^ 



(Mj/^, 0, 0, 0, ), and the sum of the four-momentum of all detected parti- 
cles in the event. Pmiss is required to be larger than 0.2 GeV/c and less 
than 0.9 GeV/c for J/ip — > Dje'^Ue and larger than 0.2 GeV/c and less 
than 1.0 GeV/c for J/ip -^ D^e^u^ to reduce background from J/ip decay 
to 4-prong final states with misidentified particles but with no missing par- 
ticles. We further require the absolute value of Umiss, which is defined as 
Umiss = Ejniss — Pmiss, to bc Icss than 0.1 GcV, to rcjcct backgrounds from 
J/ijj decaying to K^, rj, and partially missing vr" final states, which were not 
rejected by prior criteria. After all requirements, the invariant mass distribu- 



tions of D: 



K+K-TT- and D 



K' 



TX TT 



are shown in Figs. 3(a) and 



(b), respectively. 
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Fig. 3. Invariant mass distribution of (a) D^ 
and (c) W -^ K+7r-,/s:+7r-7rO,iv:+7r-7r-7r+ 
bars. The histogram is MC simulated signal events, unnormalized. 



-^ K+K't:-, (b) D- -^ K+TT'TT-, 
Data are shown by dots with error 



3.2 J/iP^D^e^ 



D^ mesons are reconstructed in three decay modes: D'^ — > K'^tt~, D^ — >■ 
i^+7r~7r°, and Z}° -^ K~^ 7r~ 7i~ Tr~^ . In D^ — >• ii'"'"7r~7r°, the 7r° candidates are 
reconstructed from a pair of photons. The 7r° decay angle \E^i — E^2\/Pit0 
is required to be < 0.9 for K^tt^tt'^ in order to remove events in which the 
TT^ is falsely reconstructed from a high-energy photon and a spurious shower. 
Here, P^o is the momentum of 7r° and E^i and E^2 are the energies of the 
two photons. Candidate events are required to have four or six charged tracks 
which satisfy the charged track selection criteria and total charge zero. Four- 



constraint (4C) energy-momentum conservation kinematic fits are performed 
under tlie J /ip -^ e~^e~K~^n~ or tlie J/ip -^ e~^e~K~^7i~n~7i~^ liypotlieses, 
and tlie xlc i^ required to be less than 15. The J/ip — ^ e"'"e~K+7r~77 hy- 
pothesis is subjected to a five-constraint (5C) fit where the invariant mass 
of the 77 pair associated with the 7r° is constrained to 171,^0, and xlc < 
15 is required. For D'^ -^ A'+tt^tt", if more than one combination of two 
good photons passes the kinematic fit, the combination with the smallest x^ 
is chosen. Gamma conversions and 7r° Dalitz decays may produce electron- 
positron pairs with low electron-positron invariant mass. The requirement 
TTiee > 0.2 GeV/c^ eliminates most of this type of background. The combined 
D^ -^ K^n^ , K^Ti^Ti^, K^Ti^Ti^ii^ invariant mass distribution is shown in 
Fig 3(c). 



4 Monte Carlo Simulation 



Monte Carlo (MC) simulation is used for the determination of mass resolutions 
and detection efficiencies. We simulate J/ip -^ D~e^Ve-, J /i' -^ -D^e+z/g, and 
J/ip -^ D^e^e~ production and decay, including the detector response. In 
these simulations, 200000 events each are generated, and the decay branching 
fractions of D~ , D~ , and D'^ are taken from the world average values [9]. The 
detection efficiencies and branching fractions obtained are listed in Table 1. 

Table 1 



ion efficiencies and branching fractions [9 




Decay Mode 


a ■ 


B^ 


D- ^(l)TT-{(P^ K+K^) 
D- -^ K-K*^{K*^ -^ K+7T-) 


5.39% 


(3.6±l.l)%x(49.1±0.6)% 


4.29% 


(3.3±0.9)%x2/3 


D- -^ K+TT-TT- 


8.67% 


(9.2±0.6)% 




5.24% 


(3.80±0.09)% 


2.31% 


(13.0±0.8)% 


2.15% 


(7.46±0.31)% 



5 Systematic Errors 



The largest systematic errors on the branching fraction measurements come 
from the uncertainties of the MDC simulation (including systematic uncer- 
tainties of the tracking efficiency) and the Umiss or kinematic fit requirements. 



The photon detection efficiency has been studied with several different meth- 
ods using J/ip -^ p°7r° decays [10]; the difference between data and Monte 
Carlo simulation is about 2% for each photon. We estimate a systematic error 
of 4% for D" — ^ K^-K^-K^ . The uncertainty for final states with no isolated 
photons is also studied using J/ip -^ pn decays and is about 2%. The system- 
atic error from electron identification is estimated to be 5% for each electron. 
The pion and kaon identification is studied, and the difference between data 
and Monte Carlo simulation is 1.5% for each charged track, which is treated 
as a systematic error. The error on the intermediate decay branching fractions 
of D~, D~, D^, 0, and K*'^ are taken from the PDG [9]. The statistical error 
of the Monte Carlo sample is taken into account. The total number of J/ip 
events is (57.7±2.7)xl0^ [11], determined from inclusive 4-prong hadronic 
final states, and the uncertainty, 4.7%, is taken as a systematic error. The 
systematic errors from all sources, as well as the total, are listed in Table 2. 

Table 2 

Summary of the systematic errors. 





J/-0 -^ Ds e+l^e 


J/^-^D-e+Ue 




J/-4> -^ D"e+e 


- 




DJ -^cpiT-, K"K*° 


D- ^ K+n-n- 


- 


DO -> K+TT 


- 


K+tt-ttO 


K+TT-n-TT+ 


MDC Simulation 


18.6% 


11.6% 


20.6% 


Gamma 


2.0% 


2.0% 


0.0% 


4.0% 


0.0% 


ePID 


5.0% 


5.0% 


10.0% 


TT, K PID 


4.5% 


4.5% 


3.0% 


3.0% 


6.0% 


BiD.s,D) 


25% 


6.5% 


2.4% 


6.2% 


4.2% 


MC Statistics 


5.1% 


2.5% 


2.7% 


number of J/4> 


4.7% 


total 


32.7% 


15.9% 


23.8% 


24.8% 


24.6% 



6 Results 



No excess of J/ip -^ Dje^Ug, J/ip -^ D^e^u^, or J /ip -^ D^e^e~ events above 
background is observed. The upper limit on the branching fractions of these 
decay modes are calculated using 



B < 



n 



obs 
UL 



Nj/,,eB(l - a-^y^) 



"J/i^ 



where n^£ is the upper limit of the observed number of events at the 90% 
confidence level. Nj/^ is the number of J/ip events, and eB is defined as 

N 

Y^EiBi, where Ei and Bi are the detection efficiency and branching fraction for 

i=l 

decay channel i and A^ is the number of decay modes, which are listed in Table 
1. For instance, for J/ip -^ D^e^Ve-, the D~ mesons are reconstructed in two 



decay modes, and A^ is equal to 2. The systematic error in the measurement 
is taken into consideration by introducing 1 — o"*^* in the denominator of the 
branching fraction calculation. 

We obtain upper limits for the observed number of events at 90% confidence 
level of 3.55 for J/V' -^ D^e+u^, 4.64 for J/ip -^ D-e+u^ and 3.07 for J/ip -^ 
D^e^e~ , using a Bayesian method with uniform prior above zero [9]. The 
likelihood distributions and the 90% CL. limits are shown in Fig. 4. The 
likelihood values for each value of the number of events are obtained by fitting 
the distributions shown in Fig. 3 with a signal shape determined from MC 
simulation and a second order polynomial to describe background. The upper 
limits are obtained from the integral of the normalized likehood at the 90% 
confidence level. The numbers used in the branching fraction calculations are 
summarized in Table 3. 

Table 3 

Numbers used in the calculation of upper limits on the branching fractions of J/ijj -^ 
Dj /D^e+Ve + c.c. and J/ip -^ We~^e~ + c.c. 





J/tp -^ Dje+Ve + c.c. 


J/V' -^ D'e+Ue + c.c. 


J/ip -^ L>Oe+e- + c.c. 


„obs 


3.55 


4.64 


3.07 


eB 


1.90x10^3 


7.98x10-3 


6.61x10-3 


Sys. Err. 


32.7% 


15.9% 


24.8% 


S(90%C.L.) 


< 4.8x10"^ 


< 1.2x10-5 


< 1.1x10-5 



In summary, we have searched for the decays J/ip -^ Dje^Ue, J/ip -^ D^e^Ue, 
and J/ip -^ D^e'^e~^ using 5.8 x 10"^ J/ip events acquired by the BESII detector 
at the BEPC e+e- collider. No evidence for any of these decays is found. The 
final results for the 90% confidence level upper limit of the branching fractions 
are given in Table 3. The upper limits on the branching fraction for decays 
J/ip -^ D~e^Ue, J /ip -^ D~e~^i'e, and J/ip -^ D^e^e~ are not inconsistent 
with the standard model. 
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